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Quasi-spherical accretion in wind-fed X-ray pulsars is discussed. At X-ray luminosities < 4 x 
10 36 erg/s, a hot convective shell is formed around the neutron star magnetosphere, and subsonic 
settling accretion regime sets in. In this regime, accretion rate onto neutron star is determined 
by the ability of plasma to enter magnetosphere via Rayleigh-Taylor instability. A gas-dynamic 
theory of settling accretion is constructed taking into account anisotropic turbulence. The angular 
momentum can be transferred through the quasi-static shell via large-scale convective motions 
initiating turbulence cascade. The angular velocity distribution in the shell is found depending on 
the turbulent viscosity prescription. Comparison with observations of long-period X-ray wind-fed 
pulsars shows that an almost iso-angular-momentum distribution is most likely realized in their 
shells. The theory explains long-term spin-down in wind- fed accreting pulsars (e.g. GX 1+4) 
and properties of short-term torque-luminosity correlations. The theory can be applied to slowly 
rotating low-luminosity X-ray pulsars and non-stationary accretion phenomena observed in some 
SFXTs. 
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1. Introduction 

Accreting magnetized neutron stars (NS) are observed as X-ray pulsars (XPSR). The observed 
characteristics of XPSRs are their X-ray luminosity L x , pulse period P* (or angular frequency 
CI* = 2n/P*) and its first time derivative P* (or ft)*), the pulse profile at different energy bands, 
and X-ray spectrum. The timing properties (i.e., the temporal behavior of ft)*) and their correlation 
with variations of L x provide important clues on the interaction mechanism of accreting plasma 
with the rotating NS magnetosphere [ffl]. Spectral features (especially, cyclotron lines) allow direct 
probing the strength of the near-surface magnetic field of NSs [||]. Their dependence on variable 
X-ray luminosity in persistent (Her X-l, [^]) and transient (e.g. 4U01 15+63, [f|]) XPSRs provides 
information on the accretion column properties (its geometry etc.). 

XPSRs can be persistent (Her X-l, Cen X-3, Vela X-l, millisecond XPSRs, etc.) or tran- 
sient (A0535+26, 4U01 15+63, etc.), and are found in both high-mass and low-mass X-ray binaries 
(HMXB: Cen X-3, Vela X-l, transients, LMXB: Her X-l, 4U1626-67, millisecond XPSRs). The 
matter from the secondary star can be supplied to NS via Roche lobe overflow (the most frequent 
case in LMXBs) or can be captured from stellar wind of the optical companion (in HMXB). In 
the first case, the specific angular momentum of accreting matter j m is usually high enought for an 
accretion disk to be formed around the NS magnetosphere, i.e. j m > Jk{Ra) = VGMRa, where 
Ra is the characteristic radius of magnetopshere (the Alfven radius), M is the NS mass, G is the 
Newton gravity constant 1 . 

2. Two regimes of quasi-spherical accretion 

In the case of accretion from the wind considered here, the accretion flow geometry can be 
more complicated: the condition for the disk formation j m > Jk(Ra) can be satisfied or not, again 
depending on the size of the magnetosphere Ra and on the specific angular momentum of grav- 
itationally captured matter from the stellar wind. The last quantity can be expressed through the 
characteristic gravitational capture radius Rb = 2GM/v 2 (the Bondi radius), where v 2 = v 2 + v 2 orh 
is the relative velocity of wind (moving with velocity v vv )and the NS (moving with orbital velocity 
v or b through the wind) as j m = kQ.i,R B . Here Q.b is the orbital angular frequency, k is the numer- 
ical coefficient which can be positive (prograde j m ) or negative (retrograde j m ) [gj], depending on 
properties of (generally inhomgeneous) stellar wind. 

There can be two different regimes of quasi-spherical accretion. The captured matter heated 
up in the bow shock at ~ Rb to high temperatures ksT ~ m p v 2 . If the characteristic cooling time 
of plasma t coo i is smaller than the free-fall time tff = Rb/ \/2GM /Rb, the gas falls supersonically 
toward the magnetosphere with the formation of a shock. This regime is usually considered in 
connection with bright XPSRs [0|. The role of X-ray photons generated near the NS surface is 
two-fold: first, they can heat up plasma in the bow-shock zone via photoionization, and second, they 
cool down the hot plasma near the magnetosphere (with k B T ~ GM /Ra) by Compton processes 
thus allowing matter to enter the magnetosphere via the Rayleigh-Taylor instability |gjp. 

'in the extreme case of very large magnetospheres, Ra ~ a, where a is the orbital radius, the flow does not form a 
disk; this case is realized in AM Her binaries (polars). 
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In the free-fall accretion regime, the X-ray luminosity (the mass accretion rate M) is deter- 
mined by the rate of gravitational capture of stellar wind at The accretion torque applied to NS 
due to plasma-magnetopshere interaction is always of the same sign as the specific angular momen- 
tum of the gravitationally captured stellar wind, i.e. the NS can spin-up or spin down depending on 
the prograde or retrograde j m . 

If the relative wind velocity v at Rb is slow (< 80 km/s), the photoionization heating of plasma 
is important, but the radiation cooling time of plasma is shorter than the free-fall time, so the 
free-fall accretion regime is realized. If the wind velocity is larger than ~ 80 — 100 km/s, the 
post-shock temperature is higher than 5 x 10 5 K (the maximum photoionization temperature for 
a photon temperature of several keV); for L x < 4 x 10 36 erg/s, the plasma radiative cooling time 
is longer than the free-fall time, so a hot quasi-spherical shell is formed above the magnetosphere 
with temperature determined by hydrostatic equilibrium [0], [|l(i|]. Accretion of matter through 
such a shell is subsonic, so no shock is formed above the magnetosphere. The accretion rate M is 
now determined by the ability of hot plasma to enter magnetosphere. This is the settling accretion 
regime. 

3. Settling accretion regime: theory 

Theory of settling accretion regime was elaborated in [jl(J. In this regime, the accreting matter 
subsonically settles down onto the rotating magnetosphere forming an extended quasi-static shell. 
This shell mediates the angular momentum transfer to/from the rotating NS magnetosphere by 
voscous stresses due to large-scale convective motions and turbulence. The settling regime of 
accretion can be realized for moderate accretion rates M < ~ 4 x 10 16 g/s. At higher accretion 
rates, a free-fall gap above the neutron star magnetosphere appears due to rapid Compton cooling, 
and accretion becomes highly non-stationary. 

Mass accretion rate through the hot shell is determined by mean velocity of matter entering 
the magnetosphere, u(Ra) = /(«) a/ 2GM/Ra- The dimensionless factor f{u) is determined by the 
Compton cooling of plasma above magnetosphere and the critical temperature for Rayleigh-Taylor 
instability to develop [||], and is found to be 

f(u)^03Mf l ^- Q l/U , (3.1) 

where M 16 = M/[10 16 g/s] and ju 30 = Ai/[10 30 G cm 3 ] is the NS magnetic moment. The definition of 
the Alfven radius in this case is different from the value used for disk accretion: 

^ io5 w(£f- (3.2) 

As the X-ray luminosity increases, Compton cooling occurs faster, and when f{u) — >xi 0.5, the 
sonic point appears in the accretion flow above the magnetosphere, and the accretion becomes 
highly non-stationary. Thus the condition f(u) = 0.5 determines the critical X-ray luminosity for 

36 1 /4 

the steady settling accretion: L x < 4 x 10 £i 3 o erg/s. 

Accretion torques applied to NS in this regime are determined not only by the specific an- 
gular momentum of captured matter j m (as is the case of the free-fall accretion), but also by the 
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possibility to transfer angular momentum to/from the rotating magnetosphere through the shell by 
large-sclae convective motions. This means that the NS spin-down becomes possible even for pro- 
grade j m . The plasma-magnetosphere interaction results in emerging of the toroidal magnetic field 
B t = K\B p ((O m — co*)tj mt , where C0 m is the angular frequency of matter at the Alfven radus, Ki the 
dimensionless coupling coefficient which is different in different sources, ti nst is the characteristic 
time of RT instability. The gas-dynamical treatment of the problem of angular momentum transfer 
through the shell by viscous turbulence stresses Jl0| ] showed that co m ps ^(Ra/Rb) 11 , where the 
index n depends on the character of turbulence in the shell. For example, in the case of isotropic 
near-sonic turbulence we obtain n ~ 3/2, i.e. quasi-Keplerian rotation distribution. In the more 
likely case of strongly anisotropic turbulence (because of strong convection) we find n^2, i.e. an 
iso-angular-momentum distribution. The NS spin evolution equation reads: 



I CO* = AM^rr _ BM 3/n , (3.3) 

where / is the NS moment of inertia. The (independent of M) coefficients A and B for the case 
n = 2 are (in CGS units): 

A ps 5.3 x 10 31 £ijui ( —Z — -) 4 ( ^t) , B ps 5.4 x lO^uM ( (3.4) 
^ 30 VlOOOkm/s/ V 10d / ^ 30 V100sy v ' 



The function ri)*(M) reaches minimum at some M cr . By differentiating Eq. (3.3) with respect to M 



and equating to zero, we find M cr 



B 3 

A (3+2n) 



"" . At M = M cr the value of CO* reaches an absolute 



minimum (see Fig. [[]). In terms of the dimensionless mass accretion rate y = -M- Eq. (3.3) can be 



written in the form 

m*=AMTy^U-{^f\ (3.5) 

where the frequency derivative vanishes at y = y^. yo = (^^) 2 " ■ The qualitative behaviour of CO* 
as a function of y is shown in Fig. ||. So in the settling accretion regime transitions from spin-up 



(y > yo) to spin-down (y < yo) occurs with changing L x . By varying Eq. ( |3.5[ ) with respect to y we 
find: 

I(8co*)=I^(8y) = ^AMTy-- (8y). (3.6) 

We see that depending on whether y > 1 or y < 1, correlated changes of 8 CO* with X-ray flux 
should have different signs. Indeed, for GX 1+4 in [JIT)], JT^ ] a positive correlation of the observed 
8P with 8M was found using Fermi data (see also jl3| ] for BATSE data). This means that there is 
a negative correlation between 8(0* and 8M, suggesting y < 1 in this source. 

4. Settling accretion regime: observations 

First, the settling accretion regime explains long-term spin-down states as observed, for exam- 



ple, in GX 1+4 [12]. Another immediate application of Eq. (B3b to observations is the estimation 



of the NS magnetic field in XPSRs with equilibrium spin period (i.e. where CO* = on average): 
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Figure 1: A schematic plot of 03* as a function of y. Torque reversal occurs at >'o. Sign of torque-luminosity 
correlations changes at y = I. 



(here we assumed n = 2 and maximum anisotropic turbulence). In wind-fed pulsars Vela X-l and 
GX 301-2, the NS magnetic field estimated in this way was found to be close to the value inferred 
from the cyclotron line measurements (see [ |Tfj| ] for more detail). 

This equation also implies that the equilibrium period of XPSRs in the settling accretion 
regime is 

, 12/11, v-4/11 ( Pb 



F «"™*®»s M * "{mjim^fs) ■ <4 - 2) 

There are known several very slowly rotating XPSRs, including some in HMXB (4U 2206+54 with 
P* = 5550 s), 2S 01 14+65 with P* = 9600 s) and some in LMXB (e.g. 3A 1954+3 19, P* = 19400 s 
[14]). Assuming disk accretion in such pulsars would require incredibly high NS magnetic fields. 
But in 4U 2206+54 the NS field is measured to be normal \^\ , B ~ 3 x 10 12 G (in 2S 0114+65 
possibly too, with B ~ 2.5 x 10 12 G [16]). X-ray luminosity in these pulsars is low, ~ 10 36 erg/s, 
which suggests the settling accretion regime in these pulsars. 

Another possible implication of the settling accretion theory can be for non-stationary phe- 
nomena in XPSRs. For example, a dynamical instability of the shell on the time scale of the order 
of the free-fall time from the magnetosphere can appear due to increased Compton cooling and 
hence increased mass accretion rate in the shell. This may even result in a complete collapse of 
the shell resulting in an X-ray outburst with duration similar to the free-fall time scale of the en- 
tire shell (~ 1000 s). Such a transient behaviour is observed in supergiant fast X-ray transients 
(SFXTs) (see e.g. [0]). Slow X-ray pulsations are found in some of them (e.g. IGRJ16418-4532, 
P* rj 1212 s [||]). In this source, the regular pulsations distinctly seen at a small X-ray luminosity 
L x ~ 10 34 erg/s were found to disappear when the average X-ray flux increased by more than an 
order of magnitude, and quasi-regular flares were observed on a time scale of 200-250 s. In such a 
low-luminosity source the magnetospheric radius should be very large, 7?^ — 5 x 10 10 cm (assum- 
ing the standard value for the NS magnetic field /I30 = 1). So the observed flaring behavior can 
be the manifestation of a Rayleigh-Taylor instability from the magnetospheric radius which takes 

3 /2 

place on the dynamical time scale ~ /? A / y GM. Note that at small X-ray luminosities the radi- 



A 

ation cooling (and not the Compton cooling) controls plasma entering the magnetosphere, so the 
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critical luminosity for unstable accretion becomes lower than 4 x 10 erg/s found for the Compton 
cooling case. 

5. Conclusion 

At X-ray luminosities < 4 x 10 36 erg/s wind-fed X-ray pulsars can be at the stage of subsonic 
settling accretion. In this regime, accretion rate onto NS is determined by the ability of plasma to 
enter magnetosphere via Rayliegh-Taylor instability. A gas-dynamic theory of settling accretion 
regime is constructed taking into account anisotropic turbulence [|1C[]. The angular momentum can 
be transferred through the quasi-static shell via large-scale convective motions initiating turbulence 
cascade. The angular velocity distribution in the shell is found depending on the turbulent viscosity 
prescription. Comparison with observations of long-period X-ray wind-fed pulsars GX 301-2 ad 
Vela X-l shows that an almost iso-angular-momentum distribution is most likely realized in their 
shells. The theory explains long-term spin-down in wind- fed accreting pulsars and properties 
of short-term torque-luminosity correlations. Long-period low-luminosity X-ray pulsars are most 
likely experiencing settling accretion too. Spectral and timing measurements can be used to further 
test this accretion regime. 
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